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We study the voter model with a finite density of zealots — voters that never change opinion. For 
equal numbers of zealots of each species, the distribution of magnetization (opinions) is Gaussian 
in the mean-field limit as well as in one and two dimensions, with a width that is proportional to 
l/\fZ, where Z is the number of zealots, independent of the total number of voters. Thus just a 
few zealots can prevent consensus or even the formation of a robust majority. 
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I. INTRODUCTION 

The voter model [l[ is one of the simplest examples of 
cooperative behavior that has been used as a paradigm 
for the dynamics of opinions in socially interacting pop- 
ulations. In the voter model, each node of a graph is 
occupied by a voter that has two opinion states, de- 
noted as + and — . Opinions evolve by: (i) picking a 
random voter; (ii) the selected voter adopts the state 
of a randomly-chosen neighbor; (iii) repeat these steps 
ad infinitum or until a finite system necessarily reaches 
consensus. Naively, one can view each voter has having 
no self confidence and thus takes on the state of one of 
its neighbors. This evolution resembles that of the Ising 
model with zero-temperature Glauber kinetics 0], but 
with one important difference: in the Ising model, each 
spin obeys the state of the local majority; in the voter 
model, a voter chooses a state with a probability that is 
proportional to the number of neighbors in that state. 

There are three basic properties of the voter model 
that characterize its evolution. The first is the exit prob- 
ability, namely, the probability that a finite system even- 
tually reaches consensus where all voters are in the + 
state, E+(po), as a function of the initial density po of + 
voters. Because the mean magnetization, defined as the 
difference in the fraction of + and — voters (averaged 
over all realizations and histories), is conserved on any 
degree-regular graph, and because the only possible final 
states of a finite system are consensus, E + (p ) = p 

A second basic property is the mean time Tjy to reach 
consensus in a finite system of N voters. For regular 
lattices in d dimensions, it is known that Tjv scales as 



N 2 



1, as AT In AT in d = 2. and as N in d > 2 



[3, [If . In contrast, Tjv generally scales sublinearly with 
A" on heterogeneous graphs with broad degree distribu- 
tions [4]. Defining pk as the k th moment of the degree 
distribution, then T/v ~ Np 2 /p2, which grows slower 
than linearly in A" for a sufficiently broad degree distri- 
bution. Finally, the 2-point correlation function G2(r), 



defined as the probability that 2 voters a distance r apart 
are in the same state, asymptotically decays as r 2 ~ d on 
a regular lattice when the spatial dimension d > 2 [1, [f| . 
This decay is the same as that of the electrostatic poten- 
tial of a point charge, a correspondence that has proven 
useful in analyzing the voter model. 

In this work, we investigate an extension of the voter 
model in which a small fraction of the population are 
zealots — individuals that never change opinion. The ef- 
fect of a single zealot § or a small number of zealots Q 
on primarily static properties of the voter model has been 
studied previously, and considerable insight has been 
gained by exploiting the previously-mentioned electro- 
static correspondence. The role of zealots has also been 
investigated in a majority rule opinion dynamics model 
Q, where again equal densities of zealots of each type 
prevent consensus from being achieved. One motivation 
for our work is the obvious fact that consensus is not the 
asymptotic outcome of repeated elections in democratic 
societies. One such example is the set of US presidential 
elections @ , where the percentage of votes for the winner 
has ranged from highs of 61.05% (Johnson over Goldwa- 
ter 1964) and 60.80% (Roosevelt over Landon 1932) to 
lows of 47.80% (Harrison minority winner over Cleveland 
1888) and 47.92% (Hayes over Tilden 1876). In this com- 
pilation, we exclude elections with substantial voting to 
candidates outside the top two (Fig. [TJ. 

This example, as well as election results from many 
democratic countries, show in an obvious way that con- 
sensus will never be achieved in large voting populations. 
This fact motivates us to to investigate an opinion dy- 
namics model in which consensus is stymied by the pres- 
ence of zealots. Because of the competing influences of 
the zealots and the tendency toward consensus by the 
voter dynamics, the magnetization fluctuates with time 
in a manner that can be made to qualitatively mimic, for 
example, the U.S. presidential election results (Fig. [IJ. 
Upon averaging over a long time period, these time- 
dependent fluctuations lead to a stationary magnctiza- 
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FIG. 1: Time dependence of the magnetization for single re- 
alizations of 1000 voters on the complete graph for Z = 2 
(black), 16 (red), and 128 zealots (blue). Data were smoothed 
over a 1% range. Also shown are U.S. presidential election re- 
sults (circles) from 1876-2004 (corresponding to t = and 
t = 3000 respectively) where the magnetization is defined as 
the difference in the vote fraction of the top 2 candidates. 



tion distribution whose properties are the main focus of 
this work. 

The basic question that we wish to address in the voter 
model with a subpopulation of zealots is: what is the na- 
ture of the global opinion as a function of the density of 
zealots? One of our main results is that equal but very 
small numbers of zealots of both types leads to a steady 
state with a narrow Gaussian magnetization distribution 
centered at zero. Here the magnetization is simply the 
difference in the fraction of voters of each species. Thus a 
small fraction of zealots is surprisingly effective in main- 
taining a steady state with only small fluctuations about 
this state. 

It should also be mentioned that there are a variety 
of simple and prototypical opinion dynamics models, in 
which lack of consensus is a basic outcome, including the 
multiple-state Axelrod model llOl. the bounded compro- 
mise model of Weisbuch ct al. [llj and its variants [l2| • 
For these models, the consensus preventing feature typ- 
ically is the absence of interaction whenever two agents 
become sufficiently incompatible. As a function of basic 
model parameters, the fraction of incompatible agents 
can grow, leading to cultural fragmentation and an at- 
tendant steady or static opinion state. 

In the next section, we define the model. Then in 
Sees. IIIII and IIV1 we solve the model in the mean-field 
limit and on a one-dimensional periodic ring. We then 
investigate the behavior on the square lattice by numeri- 
cal simulations in Sec.lVland find behavior that is quanti- 
tatively close to that in the mean-field limit. Finally, we 
conclude and point out some additional interesting fea- 
tures of the role of zealotry on the voter model in Sec. lVIl 



II. THE MODEL 

The population consists of N voters, with a fixed num- 
ber of zealots that never change opinion, while the re- 
maining voters are susceptible to opinion change. Each 
voter can be in one of two opinion states, +1 or — 1 
that we term "democrat" and "republican" , respectively. 
Thus the system consists of Z+ democrat and Z_ repub- 
lican zealots, as well as N + democrat and AT_ republican 
susceptibles. Each type of voter evolves as follows: 

1. Susceptible democrats can become republicans; 

2. Susceptible republicans can become democrats; 

3. Zealot democrats are always democrats; 

4. Zealot republicans are always republicans. 

Each agent, whether a zealot or a susceptible, has the 
same persuasion strength that we set to 1. That is, after 
a susceptible voter selects a neighbor, the voter is per- 
suaded to adopt the state of this neighbor with probabil- 
ity 1. Because the total population comprises of agents in 
one of four possible states, we have N = N + + AT_ + Z + + 
Z- . Since the number of zealots is fixed, the total number 
of susceptible individuals S = N - Z + - Z_ = N + + N_ 
is also conserved. The dynamics is a direct generalization 
of voter model and consists of the following steps: 

1. Pick a random voter, if this voter is a zealot nothing 
happens. 

2. If the selected voter is a susceptible, then pick a 
random neighbor and adopt its state; note that if 
the selected voter and the neighbor are in the same 
state, nothing happens in the update. 

3. Repeat steps 1 & 2 ad infinitum or until consensus 
is reached. 

We will investigate this model on the two natural geome- 
tries of the complete graph, a natural realization of the 
mean- field limit, and regular lattices. For the complete 
graph, all other voters in the system are nearest-neighbor 
to any voter. Thus the complete graph has no spatial 
structure, a feature that allows for a simple solution. In 
contrast, when the voters live on the sites of a regular 
lattice, a voter can be directly influenced only by its the 
nearest neighbors. 



III. DYNAMICS ON THE COMPLETE GRAPH 

On the complete graph, the state of the population 
may be characterized by the probability P(N+, N-, t) of 
finding N± susceptible voters at time t. Since N- = 
S — N + , we merely need to consider the master equation 



3 



for P(N + , t), which reads 



5=±l 



- J2 P(N+,t)W(N + ^N + + S). (1) 



5=±1 

The first term accounts for processes in which the num- 
ber of susceptible democrats after the event equals N + , 
while the second term accounts for the complementary 
loss processes where N + N + ± 1. Here W represents 
the rate at which transitions occur and is given by 



St W(N + -> N + + 1) 



St W(N+ ->N+-l) = 



N_(N+ + Z+) 

N(N - 
N+(N- 



1) 



N(N — 1) 



(2) 



The first line is the probability of choosing first a re- 
publican susceptible and then a democrat (susceptible or 
zealot), for which a susceptible republican converts to a 
susceptible democrat in the voter model interaction. We 
choose St = N^ 1 , so that, on average, each agent is se- 
lected once at each time step. 

While it is usually not possible to solve an equation 
of the form JT]), analytical progress can be achieved by 
considering a continuum N — > oo limit of the master 
equation and performing a Taylor expansion [13j . For 
this purpose, we introduce the rescaled variables n = 
N+/N, z± = Z±/N, and also s = 1 — z + — z_ so that 
s — n = N-/N. In the continuum limit, the reaction rates 
now become 

Win -> n + A -1 ) = N(s-n)(n + z+) 

W(n — > n — TV -1 ) = Nn(s-n + z-). (3) 

Expanding ([I]) to the second order in the variable n, we 
find the following Fokker-Planck equation [e| E3, EEL EH 

dP(n t) d Id 2 

dt = ^ + 2d^ lW P (n,t)]M) 

where (see e.g., Chap. VII of Ref. [l3T ]) 
a ( n ) = Sn W(n — » n + Sn) 

Sn=±l/N 

= [z+s-n(l - s)] ; 

6n=±l/N 

= [{n + z + )(s - n) + n(s + z_ - n)]/N. 

The first term on the right-hand-side of Eq. ([4]) leads 
to the deterministic mean-field rate equation h(t) = a, 
with solution 



»(*) = 



? i(0) - 



z+s 



,-(l-s)t 



(5) 



Thus an initial density of susceptible democrats in an 
infinite system exponentially relaxes to the steady-state 
value n* = z+s/(l — s). Correspondingly, the magnctiza- 
tionm = (N++Z+—N- — Z-)/N attains the steady-state 
value (z + — z_ )/ (z + + z_ ) . When the number of agents is 
finite, however, finite-size fluctuations arise from the dif- 
fusive second term on the right-hand side of Eq. (j4|) . This 
term leads to a steady-state probability distribution with 
a finite width that is centered at n* . In what follows, we 
examine these fluctuations around the mean-field steady 
state when N and Z± are both finite. 



A. Stationary Magnetization Distribution 

According to the Fokker-Planck equation (HJ), the sta- 
tionary distribution P(n) obeys 



a(n)P(n)-~\fi{n)P(n,t)] = 0, 
whose formal solution is 



P(n) 



0(n) 



(6) 



(7) 



Since the density n of agents in the state +1 ranges from 
to s, the normalization constant Z is obtained by re- 
quiring J dnP{n) = 1. This condition gives 



exP ( 2 So dn ' 
W) 



dn 



We are particularly interested in the distribution of 
the magnetization P(jn) in the continuum limit, which 
directly follows from (JT)) through a simple change of vari- 
ables. We first consider the system with the same number 
of zealots of each type, and then the asymmetric system 
with unequal numbers of zealots of each type. 



B. Symmetric case: Z+ = Z_ = Z 

When the number of zealots of each species is equal, we 
write Z+ = Z- = Z . The rate equation (|5|) then gives an 
equal steady-state density of democrats and republicans, 
n* = n + = n_ = s/2, corresponding to zero average 
magnetization, m* = 0. We now compute the stationary 
distribution of magnetization by accounting for finite-size 
fluctuations. When Z + = Z_ = Z, P(n) obeys Eq. Q 
with 

a(n) = z(l -2z- 2n), 

(3{n) = [(2n+ z)(l - 2z) - 2n 2 ]/N. 



Notice that a 



Nzd/3 
2 dn ' 



a feature that allows us to solve 



for the steady-state magnetization distribution easily. 
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FIG. 2: Steady-state magnetization distributions for 1000 
voters on the complete graph for Z — 2,8, 32, 128, and 512 
zealots (progressively steepening curves). The inset shows 
the scaled form of these distributions for Z > 8; the case 
Z = 8 slightly deviates from the rest of the distributions that 
become visually coincident. 



FIG. 3: Steady-state magnetization distributions on a com- 
plete graph of 1000 sites with unequal numbers of zealots. 
Shown left to right are the cases of (Z+,Z-) = (90,90), 
(120,60), (135,45), (144,36), (162,18). The results of voter 
model simulations and the solution to the master equations 
are coincident. The mean magnetization of the system equals 
the magnetization of the zealots: m = z + z - m 



To perform the integral in Eq. |(7J), it is helpful to trans- 
form from n to the magnetization m = (2n — s) / s which 

lies in [—1,1]. We therefore find exp ^2 J " dn' ^ = 



1 



Nz 



According to Eq. ([7]), this leads 

to the following stationary distribution of susceptible 
democrats: 



P(n) 



(zs + 2n(s — n)) 



Nz-l 



J Q S dn (zs + 2n(s — n)) 



Nz-l ■ 



(8) 



Using the fact that 2n(s — n) = ,s 2 (l — m 2 )/2, we readily 
obtain the stationary magnetization distribution: 



P(m) = -r 



z-i 



J_ 1 dm (s 1 — m 2 ) 



z-i ' 



(9) 



In the limit of large Z, we may then approximate the 
distribution by the Gaussian P(m) cx e~ m / 2<T ", with 
( t 2 = /[2 S (Z-1)]. 

When zealots are present in equal numbers, the mag- 
netization distribution quickly approaches a symmetric 
Gaussian, with a width that is inversely proportional to 
the square-root of the number of zealots and not the den- 
sity. Thus as the system size is increased, the density of 
zealots needed to keep the magnetization within a fixed 
range goes to zero. In the limiting case where there is 
one zealot of each type, the magnetization is uniformly 
distributed in [—1, 1] (Fig. [2]). Finally notice that the dis- 
tribution quickly approaches the asymptotic scaling form 
when Z > 8 (inset to Fig. [5J). 



C. Asymmetric case: Z+ ^ Z- 

Whcn the density of zealots of each type are unequal, 
we now have 

a(n) = (z + + n)s ~ n, (10) 
(3(ri) = [(2n + z + )(s-n) + nzJ\/N (11) 

in Eq. ([B]). To compute P(n) (and equivalently P(m)), it 
is now convenient to introduce the quantities 8 = z + — z- 
and r = V<5 2 + 4s. Noticing that one can write a/ (3 = 
[N(s — l)(d(3/dn) + 6(1 + s)/4]//3, one can easily com- 
pute the integral in Eq. (J7]) and thereby obtain P(n). 
Transforming from the density to the magnetization by 
n = (m + l)s/2, we obtain the following expression for 
the stationary magnetization distribution (Fig. [3]): 



ZP(m) = [1 — m(5 + ms)] 



(Z+ + Z_-2)/2 



-i (<5/2r) (2N-Z+-Z-) 



r — 8 
2 



(12) 



As in the symmetric case, Z is a normalization 
constant obtained by requiring that j_ 1 dmP(in) = 1. 
Notice that P(m) is comprised of two terms. The 
first term gives a Gaussian contribution (in the limit 
of large N) and is the analog of Eq. ([9]). The second 
term is a nontrivial contribution due to the asymmetry 
that is responsible for the skewness of P(m) which 
remains peaked around m* = z + + z ~~ ■ Close to this 
peak value, there is little asymmetry (i.e., 5 <^ 1). 
Additionally, for a large number of zealots we may 
approximate the distribution (|12[) by the Gaussian 



5 



P{m) 



\<Z- 



e -(m-m*) 2 /2a 2 



0{{m - m*)8))\, with 



IV. ONE DIMENSION 

We now turn to the one-dimensional system, where 
the behavior of the classical voter model is quite differ- 
ent from that in the mean- field limit. When zealots are 
present, however, we generically obtain a Gaussian mag- 
netization distribution, as in the mean-field case. We 
now derive the magnetization distribution — first for two 
zealots and then for an arbitrary number of zealots. 



A. Two Zealots 

Suppose that two zealots of opposite opinion are ran- 
domly placed on a periodic ring of length L. The ring is 
thus split into two independent segments of lengths Li 
and 1/2, withL = L 1 + L 2 + 2 (Fig.HJ. We take the ring to 
be large so that we may write L w L\ + L 2 . As shown in 
Fig. [U the voters in each segment coarsen and eventually 
there remains one domain of + voters that is separated 
from one domain of — voters by a single domain wall. 
Each domain wall performs a free random walk and the 
walk is reflected upon reaching the end of its segment. A 
basic fact from the theory of random walks [181 ] is that 
each domain wall is equiprobably located within the in- 
terval in the long-time limit. We now exploit this prop- 
erty to determine the magnetization distribution. 




FIG. 4: A ring divided into two independent segments by 
oppositely-oriented zealots (thick lines). Also shown is the 
state of each voter and the domain wall in each segment at 
long times (dotted lines). 

For interval lengths L\ and L 2 and respective magne- 
tizations mi and mz, the magnetization m of the entire 
ring is given by mL = m\L\ + m 2 L 2 . Thus a given value 
of to is achieved if mi and m 2 are related by (Fig. [5|) 



m 2 



mL 



L 2 



(13) 



Then the probability P(m\L 1 ,L 2 ) for a system of two 
segments with lengths L\ and L 2 to have magnetiza- 
tion equal to to is proportional to the length of the ray 



defined by Eq. (| 13[) that lies within the unit square in 
the m 1 -m 2 plane. As illustrated in Fig. [5J the distri- 
bution P<(m|Li,L 2 ), where the subscript < now signi- 
fies the range L\ < L 2 , increases linearly with m for 
— 1 < m < (Li — L 2 )/L, is constant for {L\ — L 2 )/L < 
m < (L 2 — Li)/L, and then decreases linearly with to for 
(L 2 - Li)/L < m < 1. 



m = +1 



(L 2 -L 1 



m = —1 



rn 2 



mi 





, P<(m|£i,Z 2 ) 




m 


1 (Li-L 2 ) 
1 L 


L 1 



FIG. 5: (Top) Rays of fixed magnetization (dashed) for the 
case Li < L?.. The probability for a given value of m is 
proportional to the length of the ray corresponding to this m 
value within the unit square (solid). (Bottom) The resulting 
magnetization distribution P < (m\Li, L2) for a given L\ and 
Li < L 2 . 

Using this m dependence of P < (m\Li, L 2 ) and also im- 
posing normalization, we thus find the magnetization dis- 
tribution for fixed L\, L 2 with L± < L 2 to be: 



P < (m\L 1 ,L 2 ) = { 



' L 2 {l + m) 
4Lii 2 

L 
2L 2 

L 2 (l - to) 
4LiL 2 



- 1 < m < 



\m\ < 



L 2 -L! 



< TO < 1. 



(14) 

The complementary distribution P > (m|Li,L 2 ) for L\ > 
L 2 is obtained from Eq. (|14j) by interchanging the roles 
of L\ and L 2 . 

Now we integrate over all values of L\ to find the 
configuration-averaged magnetization distribution P(m). 
The details of this calculation arc a bit tedious and are 
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P(m) 




number of + and — voters at long times is uniformly 
distributed in [—Li, Li}. We define this difference as the 
unnormalized magnetization Mj. We now make the fol- 
lowing approximations that apply when L, Z — > oo such 
that each Li is also large. In this limit, we may assume 
that each Li is independent and identically distributed. 
As a result, the sum of the unnormalized magnetizations 
over all intervals is equivalent to the displacement of a 
random walk of Z steps with each step uniformly dis- 
tributed in [—Li, Li}. 

To solve this random walk problem, we use the basic 
fact that the Fourier transform for the probability dis- 
tribution of the entire walk P(k) is simply the product 
of the Fourier transforms of the single-step distributions 
0, [TH . Since the Fourier transform of a uniform singlc- 



FIG. 6: Comparison the analytic magnetization distribution 
for two zealots on the ring (Eq. (|15)l) and simulation results 
(points). 



given in Appendix,, The final result is 



step distribution over the range [—Li 
then have 



Li} is 



sinkLj 
kLi 



WC 



r(k) = n 



sin kLi 
kLi 



(16) 



P(m) 



P < (m\L 1 ,L 2 )dL 1 + / P > (m\L 1 ,L 2 )dL 1 



l-|m|\. /1 + H 



In 



1 + |m| 



(15) 



As shown in Fig. [51 the agreement between Eq. (|15p and 
simulations is excellent. 



B. Many Zealots 

We now study the magnetization distribution when 
many zealots are randomly distributed on the ring, with 
the restriction of equal numbers of each type of zealot 
(Z+ = Z- = Z). Two distinct possibilities can arise: 

1. A segment of consecutive susceptible voters is sur- 
rounded by two zealots of the same sign. With 
voter model dynamics, these segments eventually 
align with the state of the confining zealots so that 
the segment freezes. 

2. A segment of consecutive susceptible voters is sur- 
rounded by two zealots of opposite opinion. Even- 
tually a single domain wall remains that diffuses 
freely within the segment. 

We first consider the simpler case where equal numbers 
of + and — zealots arc randomly but alternately placed 
around the ring so that no frozen segments arise. The 
segment lengths {L{\ with i = 1, 2, . . . , Z, obey the con- 
straint Y2i Li — L (ignoring the space occupied by the 
zealots themselves). 

To find the magnetization distribution, we map the 
state of the voters onto an equivalent random walk as 
follows. In a segment of length Li, the difference in the 



Since we are interested in the asymptotic limit where 
the unnormalized magnetization becomes large, we study 
the limit of V(k) for small k. Thus we expand each factor 
in V(k) in a Taylor series to first order, and then re- 
exponentiate to yield 

z 

v(k) « n^-^/e) 

i=l 

~ l-f> 2 ^/6~e- fe2 "/ 6 . 

i 

We now invert this Fourier transform to give the distri- 
bution of the unnormalized magnetization 

P( M ) = ±- [ e - k '^ L '/ 6 e- ikM dk 
2n J 



-M 2 /2a 2 M 



y/T? 



(17) 



M 



with o a M = '£ i I%/3. 

What we want, however, is the magnetization distribu- 
tion; this is related to P(M) by P(m) dm = P(M) dM. 
Wc thus find 



P(m) = 



1 



-m 2 /2al 



(18) 



where = J2i^i /3L 2 . If the number of intervals is 
large, then each Li is approximately L/Z, from which we 
obtain er^ w 1/3Z. (The result er^ = 1/3Z is exact if all 
interval lengths arc equal.) As in the mean-field limit, 
the width of the magnetization distribution is controlled 
by the number of zealots and not their concentration, so 
that a small number of zealots is effective in maintaining 
the magnetization close to zero. 

A similar approach applies in the case where the spa- 
tial ordering of the zealots is uncorrclated. In this case, 
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approximately one-half of all segments will be terminated 
by oppositely-oriented zealots and one-half by zealots of 
the same species. For the latter type of segments, the 
unnormalized magnetization will equal ±Lj equiproba- 
bly. Under the assumption that exactly one-half of the 
segments are frozen and one-half contain a single freely- 
diffusing domain wall, the analog of Eq. ([IB)) is 

z/2 . z/2 

v{k) = TT^klH J CQskL . (1Q) 

The second product accounts for frozen segments in 
which the unnormalized magnetization equals ±Lj 
cquiprobably. For these segments the Fourier transform 
of the single-step probability for a random walk whose 
steps length are ±Li equals cosfcTi. Following the same 
steps that led to Eq. (TS|) . we again obtain a Gaus- 
sian magnetization distribution, but with <r^ given by 
^ = E i 2i?/3L 2 ^2/3^. 



V. TWO DIMENSIONS 

In the classical voter model, the two-dimensional sys- 
tem is at the critical dimension so that its behavior de- 
viates from that of the mean-field system by logarithmic 
corrections. In the presence of zealots, however, the be- 
haviors in two dimensions and in mean field are quite 
close, as illustrated in Fig. [7] 



P(m) 




m 



FIG. 7: Comparison of simulations for the magnetization dis- 
tribution in two dimensions (dashed) with the mean-field re- 
sults (solid curves). The simulations are for 1000 voters with 
2, 8, 32, 128 and 512 total zealots, with equal numbers of each 
type. 

Our results for two dimensions arc based on numerical 
simulations. In our simulations, we pick a random voter 
and apply the update rules of Sec. [TTJ The unit of time 
is defined so that a time increment dt — 1 corresponds to 
N update events, so that each voter is updated once on 
average. The system is initialized with each voter equally 



likely to be in the + or the — states. From the N + vot- 
ers in the + state, Z + of them are designated as zealots, 
and similarly for voters in the — state. After the system 
reaches the steady state, we measure steady-state prop- 
erties at time intervals AT. The delay time T to reach 
the steady state depends on the lattice dimension and 
the zealot density, while AT is the correlation time for 
the system in the steady state. By making measurements 
every AT steps, we obtain data for effectively uncorre- 
cted systems. Typically, for a given initial condition, we 
made 100 measurements and then averaged over many 
configurations of zealots. 

The resulting data for the magnetization distribution 
is typically noisy, and we employ a Gaussian averaging of 
nearby points to smooth the data. If m, denotes the i th 
magnetization value, then the smoothed magnetization 
distribution at m, is defined as 

d 

P(m t )= £ P(m l+k )e-^ 2 , 

k=-d 

where the sum includes the initial point, as well as the 
d points to the left and to the right of the initial point, 
with d typically in the range 20-40. Such a smoothed 
distribution is the quantity that is actually plotted in 
Figs. [H [31 [3 and in the spatial averaged distribution in 
Fig. E 

VI. DISCUSSION 

We have shown that a small number of zealots in a pop- 
ulation of voters is quite effective in maintaining a steady 
state in which consensus is never achieved. When there 
are equal numbers of zealots of each type, the steady- 
state fraction of democrats and republicans equals 1/2; 
cquivalently, the magnetization equals zero. For unequal 
densities of the two types of zealots, the steady-state 
magnetization is simply m* = (Z + — Z-)/(Z+ + Z-), 
where Z+ and Z- are the number of zealots of each type. 
The magnetization distribution is generically Gaussian, 
P(m) oc e -{m-m') 2 /^\ w ith a cx 1/^/Z, and Z = 
Z + + Z- is the total number of zealots. A Gaussian 
magnetization distribution arises universally in one di- 
mension, on the square lattice (two dimensions), and on 
the complete graph (mean-field limit). One basic conse- 
quence of this distribution is that as the total number of 
voters N increases, the fraction of zealots needed to keep 
the magnetization less than a specified level vanishes as 

i/Vn. 

There are several additional aspects of the influence 
that zealots have on the voter model that are worth point- 
ing out. Although the time to reach consensus is infinite 
because this state can never be achieved, one can ask 
for the time until a specified plurality is first achieved. 
Equivalcntly, we can ask for the probability that the mag- 
netization first reaches a value m, when the system is ini- 
tialized with m = toq. From the above generic Gaussian 
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form of the magnetization distribution, we expect that 
the mean time for a symmetric system to first reach a 
magnetization m will thus scale as e am z , where a is a 
constant of order one. Thus one must wait an extremely 
long time before the system achieves even a modest devi- 
ation away from the zero-magnetization state when the 
number of zealots becomes appreciable. Perhaps this 
trivial fact is the underlying reason why so many demo- 
cratic countries are characterized by small majorities in 
governance. 
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FIG. 8: Comparison of simulations for the magnetization dis- 
tribution in two dimensions when the two zealots are adjacent 
(curve with peaks near ±1), maximally separated (dots), and 
averaged over many different zealot configurations (dashed). 

Another interesting feature is the role of the zealots' 
spatial positions on the steady state. For example, if 
there are only two zealots that are adjacent, one might 
expect that the effect of this "dipole" would be weaker 
than that of two separated monopolcs. This is pre- 
cisely the effect that is is observed in Fig. [5J When the 
two zealots are adjacent, their effects arc substantially 
screened and the magnetization distribution is peaked 



near m = ±1. That is, the voters show a preference for 
consensus in spite of the zealots. On the other hand, 
when the zealots are maximally separated, the magneti- 
zation distribution is close to the distribution that arises 
when averaging over possible positions of the two zealots. 

Zealots are also quite effective in reducing the total 
number of opinion changes in the system. If the popula- 
tion is close to zero magnetization, each voter typically 
has equal numbers of neighbors of each type. If the voters 
are not strongly correlated, each voter would change its 
state at a rate that is approximately equal to 1/2. How- 
ever, simulations on the square lattice show that the flip 
rate of each susceptible voter is considerably smaller. For 
example, for 1000 voters with 10 zealots (5 of each type), 
the rate of opinion changes of the susccptiblcs is around 
1/5 and this rate decreases as the density of zealots de- 
creases. 

Finally, a slight embellishment of our model could ap- 
ply to real voting patterns in a democracy with strong 
regional differences. Here it is natural to partition a pop- 
ulation into enclaves, with an imbalance of one type of 
zealot over the other in each enclave. Such a spatial dis- 
tribution would correspond to red (republican) and blue 
(democrat) states in the parlance of US electoral politics. 
It would be interesting to study if such an extension can 
actually account for election results. 

Added Note: As this manuscript was being written, we 
became aware of a very recent cprint by Chinellato et 
al. [Hj]; they study essentially the same model as in this 
work, but with a somewhat different focus than ours. 
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APPENDIX: MAGNETIZATION DISTRIBUTION FOR TWO ZEALOTS 



We want to compute the integral 



P(m) = - 

Li 



— L 

Jj P < (m\L 1 ,L 2 )dL 1 + P > {m\L l ,L 2 )dL 1 



(A.l) 



Since P < (m\Li, L 2 ) and P > (m\Li, L 2 ) have different forms in different parts of the interval [—1, 1], each of the above 
integrals needs to be split into two parts. For P<(to|I/i, L 2 ) and assuming that m > 0, the linear ramp part of the 
probability distribution needs to be used for [L 2 — L±)/L < m, which translates for L x > L(l — m)/2. Similarly, for 
P > (m\Li, L 2 ) and again for to > 0, the linear ramp must be used when {L\ — L 2 )/L < to, or L\ < L{1 + m)/2. Thus 
the above integral becomes 



P(m) 



1 
L 



r(x- 



LdL x 
2(L-Li) 



(1 - m)L 2 dLi 
4Li(L-Li) 



(1 - m)L 2 dLi 
±L X (L-L X ) 



(X+m) 



LdL x 
2Li 



(A.2) 
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Each of these integrals is then elementary. We also obtain the result for to < by reflecting the result of the above 
integral about to = to give Eq. (fT5)) . 
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